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Abstract 

The connection between continued fractions and orthogonality which is famihar for 
J-fractions and T-fractions is extended to what we call /^-fractions of type I and II. 
K^ I These continued fractions are associated with recurrence relations that correspond to 

multipoint rational interpolants. A Favard type theorem is proved for each type. We 
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then study explicit models which lead to biorthogonal rational functions. 
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1. Introduction. There is a close connection between the theory of orthogonal polyno- 
mials and continued fractions. This connection is clearest in the case of the classical moment 
problem where orthogonality is with respect to a positive measure whose Stieltjes transform 
is represented by a positive definite J- fraction (see [^, |^). This connection is known to 
extend to the quasi-definite case where the measure and the J-fraction are no longer positive 
and positive definite, respectively [0, |T^|. Another extension is by means of T-fractions and 



their connection with the trigonometric moment problem |^, |^T|, [^ . With both types of 
extensions one has a Favard type theorem establishing orthogonality with respect to a unique 
moment functional, given the three term recurrence relation satisfied by the polynomials under 
consideration. More importantly there is always a spectral measure(s) associated with these 
problems and by determining the regions where the associated continued fraction diverges 
one can locate the support of the corresponding spectral measure(s). Furthermore by finding 
the limits to which the associated continued fraction converges one can compute, at least in 
theory, the spectral measure(s). 

In this paper we shall examine two further natural extensions which, from a different 
point of view, are associated with multipoint rational interpolants or Fade approximants 
(lig, 0, [gg, ig, |2§). We call the new types of continued fractions "i?-fractions". We 



establish a Favard type theorem for /^-fractions. We also establish the existence of a natural 
Borel measure associated with i?-fractions. In general we show the existence of rational 
biorthogonality with respect to this Borel measure. We demonstrate by means of examples 
how our general approach can yield explicit biorthogonal systems of rational functions and 
can lead to the evaluation of some new types of beta integrals. 

In Section 2 we consider what we call continued fractions of type Rj. Such continued 
fractions are associated with the polynomial recurrence relation 

(1.1) P„(x) - (X- Cn)Pn-l{x) + A„(x - a„)P„_2(x) = 0. 

In addition to proving a Favard type theorem for i?7-fractions in Section 2, we also establish 
a biorthogonality relation for rational functions. Section 2 contains, as an application of our 
results, a new evaluation of a g-beta integral of Ramanujan and an alternate derivation of a 
biorthogonality relation due to Pastro |2^. Pastro's original proof involves the aforementioned 



integral of Ramanujan. It is worth noting that we obtain both the evaluation of the g-beta 
integral and the biorthogonality relation from the same result on i?/-fractions. An additional 
result is the evaluation of the Herglotz transform of the Borel measure involved. 

In Section 3 we consider continued fractions of type Ru which are associated with poly- 
nomial recurrence relation 

(1.2) Pn{x) - (X- Cn)Pn-l{x) + A„ (x - a„)(2; - bn)Pn~2{x) = 0. 



In Section 3 we follow the same plan as in Section 2. We prove a representation theorem 
for moment functionals associated with (1.2) and study the convergence of the corresponding 
continued fraction. We also identify biorthogonal rational functions that arise naturally from 
the recurrence relation (1.2). 

The letter R in our terminology is used to emphasize that the corresponding continued 
fractions are associated with rational interpolation and rational biorthogonality. 

To illustrate our ideas we include seven examples where we have applied our general 
results to concrete situations. We have already mentioned that in Section 2, we rederive a 



result of Pastro |26| . In Section 3 we give two new systems of biorthogonal rational functions 
associated with continued fraction of type Rn. The first is based on the Rn recurrence relation 
with constant coefficients which corresponds to Chebyshev polynomials. The second uses 2-F1 
contiguous relations and yields a weight function which is the integrand in the Cauchy beta 
integral. Section 4 contains two examples of biorthogonal 4^3 rational functions based on 
solutions to the Askey- Wilson three term recurrence relation. The first rederives a system of 
biorthogonal rational functions due to Al-Salam and Ismail ^. The second turns out to yield 



special cases of a biorthogonality related to g~^-Hermite polynomials, [0. Section 5 contains 
two additional examples based on a modification of the Askey- Wilson recurrence. The first 
is new and unexpected. It provides an integral representation for a 20i basic hypergeometric 
function and a special case of it leads to rational biorthogonality on [—1, 1] based on the 
elementary integral 

2 /•! VI - a;2 dx 



Ti J -1 {I - 2ax + a^){l - 25x + d"^) I - a5' 

The second, which is originally due to M. Rahman [^, evaluates a new g-beta integral of 
the type studied by Askey and Wilson and derives an associated system of rational func- 
tions biorthogonal on [—1,1]. The biorthogonality is between a {403} and {e^s} pair. The 
methodology followed in the examples combines contiguous relations for hypergeometric or 
basic hypergeometric functions with Pincherle's theorem [^], p2[. This methodology proved 



fruitful in other types of continued fractions, [Tj, |]T2|, |T6||, [|j, |^, |33 



We shall mostly follow the notation of Gasper and Rahman in Q. The shifted factorials 
and multishifted factorials are 

n 

(1.4) (a)o:=l, (a)„:=n(« + J-l), ^ > 0, 
and 

k 

(1.5) (ai,a2,- ■■,afc)n := n('^i)«' 



respectively, and a hypergeometric function is 



;i.6) rFs{ai, 02, ■ ■ ■ , ttr] h, 62, " " " , bs] Z) = rFs 



O'li Q'2, • • • , Ctr- 

bi,b2,---,bs 



E 



(Oi, a2, ■ ■ ■ , Oj-jn ^'^ 



The g-shifted and multishifted factorials are similarly defined by 



(1.7) (a;g)o:=l, (a; g)„ := ]J(1 - ag^ ^), n>0,orn = oo, 



1.8) (ai,a2,---,afc;g)„ := n('^i5 9)n- 



A basic hypergeometric series is 



(1.9) r0s(ai, 02, ■ ■ ■ , «r; &1, h, ' ' ' ', Q, z) = r<l)s 



bi,b2,---,bs 



q,z 



(ai, a2, ■ ■ ■ , a,; g)„ [(_i)"g"(„-i)/2]i+s^r^. 



:=y 

^0 {bi,b2,---,bs]q) 



(g;g)r 



Another useful notation is the short hand notation for a very well-poised basic series. We 
shall use 



(1.10) 8W7{a;b,c,d,e, f;z) := 8(1)7 



a.^faq, -^/aq, b, c, d, e, / 
a, — v^, aq/b, aq/c, aq/d, aq/e, aq/f 



q,z 



The biorthogonality based on (1.3) which we alluded to earlier is 



2 /•! 



(a(5)"(g,g;g)„ 



(1-11) -/ fm{x;S,a)fn{x;a,6)Vl-x^dx= r 'T\'''/r' ""t 2v\ '^"^."' 
with maa;{|a|, \6\} < 1 and /„ is defined by 
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(1.12) /n(^;«,^)-Y3^^^-^403 






2ax + q;2 \ q;(5, gae*'', qae 



g, g , X = cos 6*. 



The functions in (1.12) can be expressed in the simpler form 

(1.13) Ux;a,5) f ^^'^^" Mg";?)^ {-1)^^'''^' 

t^Q{q]q)k{q]q)n-k {aS;q)k I - 2aq^x + a^q'^'^' 

The expression , . -. (' n — is the Gaussian binomial coefficient. The functions {/(x; a, 5)} are 
orthogonal with respect to the weight function of the Chebyshev polynomials of the second 
kind and may play an important role in a future theory of biorthogonal rational functions. 

2.i?/-fractions. We begin this section with a Favard type theorem for the polynomial 
recurrence (2.1) below. We then outline how orthogonality may be realized in terms of the 
properties of what we will call an Ri type continued fraction. We also include an example to 
illustrate this approach. 

We first establish a Favard type theorem which is Theorem 2.1. Consider a system of 
monic orthogonal polynomials generated by 

(2.1) Pn{x) = {x- Cn)Pn-i{x) - A„ (o; - a„) P„_2(a;), 



P_i(x) := 0, Po{x) := 1, 

where 

A„+i 7^ 0, P„(a„+i) 7^ 0,n = 1,2,.... 
The recurrence relation in (2.1) can be renormalized to yield the rational recurrence relation 

(2.2) {X - a„+i) Rn{x) -{X- Cn)Rn-l{x) + \nRn-2{x) = 

with the same initial conditions, namely 

(2.3) Ro{x) := I, R^i{x) := 0. 

The renormalization is given explicitly by 

n 

(2.4) Rnix) := Pnix)/Y[{x - afc+i). 

k=l 



Theorem 2.1 Associated with the recurrence relation (2.1) there is a linear functional C 
defined on the span of {a;'^-R„(x)}^^^Q mapping it into C, and normalized by C[l] = Ai 7^ 0, 
such that the orthogonality relation 

(2.5) C[x^Rn{x)] = 0, < A; < n, 

holds. Furthermore the functional values £[a;"] and C\\Xk=i{x — cik+i)"^], n = 1,2,... are 
uniquely determined in terms of the sequences {a„+i, c„, A„ : n = 1, 2, ■ • ■}. 

Proof. Define a linear functional C whose action on Rn{x) and x"'~^Rn{x) is given by 

(2.6) C[l] := Ai, and C[Rnix)] = £[x"~^i?„(x)] = 0, for n > 1. 
We first establish (2.5) for 1 < A; < n — 1. Clearly (2.6) yields 

(2.7) = C[Ri{x)] = C[xR2ix)] = Clx'^R^ix)]. 

Apply C to (2.2) with n = 3 to obtain C[xRs{x)] = 0. Thus (2.5) holds for 1 < n < 3. When 
n > 3, formula (2.5) follows from (2.2) and (2.6) by induction over n. To prove the latter 
part of theorem observe that C[Pi{x) / {x — 02)] = implies C[l + Pi{a2)/{x — 02)] = 0. Using 
C[l] = Ai we then obtain 

C[{x-a2)-'] = -Xi/Pi{a2). 
Similarly from C[P2{x)/{{x — a2){x — 03)}] = we find 

C[l/{{x - a2)ix - 03)}] = Ai(A2 + 03 - a2)/{P2ia3)Piia2)}. 
We next use £[x^-R„(x)] for n = 1, 2 and k < n to obtain 

C[P2{x)/{x - 02)] = 'C[-A2 + {x - C2)Pi{x)/{x - 02)] = C[x - Ci - A2] = 0. 

Thus C[x\ = (ci + A2) Ai. We now continue by induction on n. Thus for each new n we use 
C[Rn{x)] = and we obtain a value for i3[l/ 11^=1(1 — cik+i)] while from C[x^~^Rn{x)] = 
and the recurrence relation (2.2) we evaluate £[a;"^^]. This establishes our theorem. 

Corollary 2.2 We have 

(2.8) C[Pn{x) Rn{x)] = £Ki?„(x)] = Ai A2 . . . A„+i, n > 0. 



Proof. Multiply (2.2) by x"~^, apply L and then use the orthogonality relation (2.5). This 
yields the two term recurrence relation £[a;"~^i?„_i(x)] = A„£[a:;"'~^i?„_2]. Taking into account 
the initial condition C\\\ = Ai we establish (2.8) and the proof of the corollary is complete. 

Note that when a„ = for n > 2 we have a Favard theorem associated with general 



T-fractions with all the moments /^[x""], n = 0, ±1, ±2, . . . uniquely determined ||T5[ ||20 |. 

Our next result is a representation theorem for the functional C based on the properties 
of a continued fraction represented as an integral transform. Recall that the polynomials of 
the second kind associated with the recurrence relation (2.1) are given by 

(2.9) Qnix) - {x- Cn)Qn-iix) + A„(x - a„)Q„_2(x) =0,n>2, 

with 

Qo{x) := 0, Qi{x) := 1 

The ratio Qn{z)/Pn{z) is the nth convergent of the infinite continued fraction 

1 \2iz- aa) A3(z - 03) 



:2.10) Ri{z) 



Z — C\ — Z — C2 — Z — C3 — 



which we assume becomes a finite fraction when z = ak,k > 2. This leads to the following 
definition. 

Definition. A continued fraction of the type (2.10) will be referred to as an /^/-fraction 
provided that it terminates when 2; = a„, ra > 2. 

We will henceforth assume that Ri{z) converges to a function which vanishes at infinity 
and whose singularity structure is given by a finite number of branch cuts and a denumerable 
number of poles so that 

da(t) 



:2.11) Rj{z):=f^, zeC\V. 

JT Z — t 



Here the measure da{t) and the multiple contour F are taken in the following generalized 
sense. In order to accomodate a pole at Zk with multiplicity m^ and residue Rk we would 
include in the right side of (2.11) a term 

Ru dt 



2TTi J\t^Zk\>e, \z-t\=e {z - t) {t - Zk) 



m,fc 



with a suitably chosen positive e. For a branch cut along a contour Tj with discontinuity a' (t) 
we would include a term 

2m Jv, z-t ' ^ ^' 
We will also assume that the domain 

(2.12) V={ul,Ti)U{zk:k = l,2,...}-, 

where {■}^ denotes a closed set. In other words we assume that the continued fraction 
converges except at the singular points of the function that it represents. 



From Pincherle's theorem [21] we then have the existence of a special minimal solution 



X^™*")(z) satisfying the same recurrence relation as Pn{z), namely 

(2.13) X(™")(^) -{z- c„) xiTH^) + A4^ - an) 4':T^(^) = 0, z^V 

but with the minimality condition 

lim Xi"''''\z)/PJz) =0, z4V. 

Hence the additional representation 

Xlr'''\z) r dajt) 

Ai(^-ai)xiT")(z) ~ Jr z-t 



:2-i4) Ri{z)= . ^ ' ..\Ln),. = L^^^ zeC\V. 



We will also assume the normalized asymptotics 

-| oo 

(2.15) Ri{z) = - +Y.dnZ-'', kl^oo, 

^ n=2 

SO that j°°^da{t) = 1. 

A further technical assumption will be that 

(2.16) {a„ : 2 <n}nP = 0. 

From the viewpoint of multiple point rational interpolants this means that the sequence of 
interpolation points associated with Ri{z) is {oo, 02, oo, 03, oo, . . .} with the a^'s distinct from 



the singular points of V, [28 



Remark. Note that Ri{z) in (2.10) does not depend on Ai,ai while they occur in (2.14). 
This seeming dependence on Ai,ai can be interpreted in two ways. Firstly the denominator 
\i{z — ai)Xl™" (-z) is determined through the recurrence relation (2.13). That is 

Ai(^ - ai)xi7")(^) := {z - Ci)X^™*")(^) - xi™")(^). 

Secondly, in explicit models where the sequences {A„}, {on} and {X^™*")(2;)} have an explicit 
analytic dependence on n, there are natural choices for Ai,ai and X^l [z). It often turns 
out that, for these natural values, Ai = 0. In these cases (2; — a„)X„™2 (-2) is singular at n = 1 
and the product \i{z — ai)X_l {z) is indeterminant. For these cases it is convenient and 
correct to define \i{z — ai)X_7 (-2) through 

\i{z - aOXl^*"^^) := lini A„(z - a„)X;™2"''(^). 

This is also true for (2.18). A similar situation occurs for i?//-type fractions in (3.16) and 
(3.17). 

We now state and prove a representation theorem. 

Theorem 2.3 Consider the three term recurrence relation (2.1) and assume the representa- 
tion (2.11) together with conditions (2.12), (2.15) and (2.16). Then the linear functional of 
Theorem 2.1 has the representation 



(2.17) C[f] = J^f{t)da{t). 



Proof. Let Ai = 1. Then /p(iQ;(t) = 1 = Ai. It remains to prove that Jj^t'' Rn{t)da{t) = 
0, < k < n. This follows from the Lemma below by taking the z ^ 00 asymptotics which 
yields for < k < n 

z-^ f t" Rn{t) da{t) = z'-^'-^Xi A2 ■ ■ ■ A„+i [1 + 0(1/ z)], 



where we have used the fact that X^^^\z)/[Xi Xi7™ (-2)] ~ IVj=2 -^j as z ^ 00 which follows 
from (2.13)-(2.15). 

Lemma 2.4 . Consider the three term recurrence relation (2.1) and assume that the repre- 
sentation (2.11) together with the conditions (2.12), (2.15) and (2.16) hold. Then we have 

A,,YiT"'(2)n"i'(2-ai) ^r z-t ' - - 



Proof. The singularities in (2.18) come from the zeros of Xi{z — ai)X_l^^ (z) and the 
discontinuities in X^^'^\z) / \i(z — ai)X_™" (z) across the contours Tj, that is, the same 
singularities as in (2.14). The singularities which would seem to appear from the denom- 
inator factor 11^=1 (-2 ~ (^j) are, without loss of generality, canceled by corresponding zeros 
of X^^^\z) (or alternatively poles of AiXi™" (-z) if one chooses a different normalization). 
In order to justify the right side of (2.18) it remains to compute the "weight" of the sin- 
gularity. If AiXi7"^(zj) = 0,xi™")(zj) ^ then it follows from (2.12) and (2.1) that 
X^^'^\zj) = Pn{zj) Xq"'\zj) since each side of this equality satisfies the same recurrence 
relation and the same initial conditions. Thus for a pole at Zj of multiplicity rrij in (2.18) one 
has a residue 



n+l / 

^,!4P„(^,)X^"H^.)/[Al^(^/-«0(^; 



m,!4P„(z,)xr"^(^,)/[Ai n(-/ -«o ( ^^^rM (^,)] 



as compared with the residue 

in (2.14). Hence we see that the pole singularities in (2.18) have a residue with an additional 
factor ZjRn{zj) as required. The situation is similar for the absolutely continuous contribution. 
Thus ii A{X^^"'\z) / Xi{z — ai)X_l (z)) is the discontinuity across a contour passing through 
the point z then A(X^™")(z)/Ai(;2-ai)xi7")(2)) = Pn{z)A{xt'"\z)/Xi{z~ai)X^^^'''\z)), 
since each side of this equality again satisfies the same three term recurrence and the same 
initial conditions at ra = —1, 0. The condition < k < n is required so that (2.18) has at least 
0(1/ z) asymptotics and thus a zero at infinity. 

We now illustrate the above theory with an example which is the g-analog of Jacobi- 
Laurent polynomials JTSf . Although our example involves a special type of i?/-fraction, namely 



a general T- fraction, we consider it to be of intrinsic interest. 

Example 2.1: 20i Functions. From the contiguous relations for 20i hypergeometric 
functions in [jl^ , or by expanding and equating powers of z, it can be verified that 



(2.19) X„+i(.)-l. + g ^^_^qn^i^)Xn{z) + q ^(i _ ^^n+l)(l _ ^^n)^n-i(^) - 



has solutions 



(2.20) A-W(.-) := ^" |°„C Ift^l:^}: ,MlM'-':k"^';,,azg^") 
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and 



(2.21) Af)(.) := ^"" '(°!ll,X"l''|: .M,IK,''-\a,''-^^,,M'V,-^ 



and a polynomial solution 

(2.22) P„(z)-^"^'^^'^^" 



,0i(g-",ag;g^-76;g,;.gVV6). 



(«?; (l)n 
The large n asymptotics of these solutions is easily seen to be given by 

(2.23) X(i)(z) ^ z- {zq'l^- q)J{azq^l^- q)^, 

(2.24) Xi'\z) ^ q-'\q'/'/z- q)oo/{hq'/' / z- qU 



and 



(2.25) P„(z) 



n/2 {b,azq^/^;q)a^ lrl^lnlV2 



II I n /9 



,n (hq^'^/z,qb;q)a 

Thus the minimal solution to the recurrence relation (2.19) is given by 



(2.26) X^™")(2) 



X«(z), |z|<|g|V2 



Pincherle's theorem |]21| then establishes the continued fraction representation 
' g-^/2^20i(l/a,g;g&;g,^ag^/^), \z\ < W^ 
z'^ 20i(g/&, g; ag^; g, hq^i'^jz), \z\ > jgl^/^ 



(2.27) i?,(^) 



where 



(2.28) Rj{z) := [z - Ci - K^^iXk z/{z - Cfc)}]"' , 



with 

(2.29) c„ = -g 



u A^-bq--'] 



1 - ag" 



A„ = q'/ 



i/2(i-g"~')(i-a&g'^"') 



(1 -ag")(l -ag'^-i) 
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Recall the Heine transformation p, (III.l)] 

(2.30) 201 (a, b; c; q, z) = — ^ — ^—^ 2<Pi{c/b, z; az; q, b). 

[c, Z] q)co 

From the above transformation we see that the 20i's on the right side of (2.27) have singu- 
larities at z = g~"~^/^/a and z = feg""*"^/^, n = 0, 1, ■ ■ ■. However, because of the respective 
conditions \z\ < \q\^^'^, \z\ > Iq]^^"^ the right side of (2.27) has no pole singularities if \aq\ < 1 
and |6| < 1. Thus when \aq\ < 1 and \b\ < 1 we also have the representation 



(2.31) R,{z)=f ^^^^, \zWil- 



with a' given by 

(2 32) olit)-= hm 1 / x(-)(t,) _ X(-)(t_) \ ^ ^ , , ^^ 

^-o+,n^o27rzA„+ii lxiT^(t+) xiT^(t_)r 



t±:=(|g|i/2±e)e*^0<^<27r. 
A calculation using (2. 26), (2. 30) and the g-Vandermonde nonterminating sum [y, 11.23] yields 
i {q^/H,q^l'^/t,q,abq]q)oo 



^2.33) a'{t) 



27rgi/2 (ogi/2t^ bq^^/t, aq^ , b; q)^ ' 



From our general theory associated with Theorem 2.3 or the general theory of T-fractions pO 
we may then state the following results 

(2.34) / r^Pn^it) a'(t) dt = 0, < k < m, 

7|t| = |g|l/2 

(2.35) / P,Ma-{t)dt=^-p^>4^. 

y|t|=|g|V2 {aq,aq^;q)m 

and 

(1 -aq)q''^/'^{abq,q;q)m 



We now derive a rational biorthogonality using (2.34) and (2.36). As a first step we shift the 
contour of integration in (2.34) and (2.36) to the unit circle |t| = 1 in order to obtain 

(2.37) f'\-''P^{t)f{a,b,t)de=^^^^^^6k,„., 0<k<m, t = e'' 
Jo [aq, bq; q)m 
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with /(a, b, t) given by 

1 iq^/H,qyyt,q,qab;q), 



f{a,b,t) 



2tt {aq^/H,bq^/'^/t,qa,qb]q)c 



The condition for shifting the contour is that no poles of a'{t) are crossed as one goes 
from |t| = \q\^^'^ to \t\ = 1. Thus one needs |g~"'+-'^''^| > \a\ and \bq^\ < 1, n = 1,2,---, 
i.e. \a\ < Iq]^^^"^ and |6| < 1. The latter two conditions are satisfied since we assumed that 
\a\ < \q\^^ and |6| < 1. 

We now take the complex conjugate of (2.37) followed by the replacements (g, a, b) — >■ 
(g, b, a) to obtain 

(2.38) r t' QUm f{a, b, t) de = i^' ""[^^ '^\"' 6k,m, < k < m, 
Jo [aq, bq; q)m 

where Qm{,z) is the polynomial 

(2.39) Q.^{z) = fp^}^ ,^^(q-m^ ^q. ^i-/^. ^, ,qi/2/^y 

[bq; q)m 

From (2.37) and (2.38) we finally obtain the biorthogonality relation 

(2.40) r Pm{t) Qn{l/t) f{a, b, t) de = i^'^'^^'^l'" 5^,„, \a\ < 1, |6| < 1. 
Jo [aq, bq; q)m 

This biorthogonality had been previously obtained by Pastro [^ using other methods. In 
[|I| the biorthogonality (2.40) was used to derive a biorthogonality for 4^3 rational functions. 
Here we have shown that the biorthogonality (2.40) is a byproduct of the theory of Rj fractions 
and we have derived a more general transform given by (2.18); namely if \aq\ < 1, \b\ < 1, 
< k < n then 

(2.41) -/ t^-"Pn(t) iqy%q'/yt;qU ^^ 



27r 7|ti=|g|i/2 z-t {aq'^/H,bq^^'^/t;q)c 

^'{?S;SS^20i(l/a,g'^+^6g"+^g,aV/'), \z\ < |g|V2 

where Pn{z) is as in (2.22). 

13 



The case z = k = n = Om (2.41) is equivalent to the original g-beta integral of Ramanujan 
with which Pastro started [^ . 

3.i?//-Practions. The format of Section 2 is repeated with first a Favard type theorem, 
then its realization in terms of the properties of an Rjj type continued fraction. We include 
an example of a new system of biorthogonal rational functions. In Section 4 we cast the recent 
results of ||I[] and |T^ in the language of Rjj fractions and use the theory outlined in this work 
to give new derivations of these results. We also discover a new set of biorthogonal rational 
functions in Section 5. 

Consider the system of polynomials {P„(x)} satisfying the recurrence relation 

(3.1) Pnix) - {x-Cn)Pn^lix) + \n(x -an){x -bn)Pn-2{x) =0, U > 1, 

and the additional assumptions 

P_l(x) = 0, Po{x) = 1, A„+i ^ 0, Pn{an+l) y^ 0, P„(6n+l) ^0, U > 0. 

The recurrence relation (3.1) can be renormalized to yield the rational function recurrence 
relation 

(3.2) {X - Qn+l) {X - bn+l) Sn{x) - {x - Cn)Sn-l{x) + \n Sn-2{x) = 0, n > 1, 

with 

n 

(3.3) Sn{x) = Pn{x)/ n [(^ - ak+i){x - hk+i). 

k=l 

We now come to an analog of Favard's theorem. 

Theorem 3.5 Given the recursion (3.1) there is a linear functional C defined on the span 
of the rational functions {x^Sn{x) : Q < k < n < oo}, mapping it into C and normalized by 
C[l] = Nq, C[xSi{x)] = Ni such that the orthogonality relation Clx'' Sn{x)] = 0,0 < k < n 
holds. Furthermore the values of >C[n"=i(a; — %+i)^^ nfcLi(a; — ^j+i)^^] for m,n = 0, 1, ■ ■ ■, 
are uniquely determined. 

Proof. Define a linear functional £ by requiring it to satisfy £[1] = A^o, C[xSi{x)] = Ni, 
C[Si{x)] = and C[Sn{x)] = C[xSn{x)] = 0, n > 2. The recurrence relation (3.2) then 
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yields C[x''Sn{x)\ = 0, 1 < A; < n. From C{xSi{x)\ = Ni and £[51 (x)] = we obtain 
C[{x — Ci)/{x — 02)] = '^[{x — Ci)/{x — 62)] = ^1 and hence 

C[l + (02 - c,)/{x - 02)] = £[1 + (&2 - c,)/{x - 62)] = N,. 

With C[l] = No this iniphes C[{x - 02)^^] = {Ni - No)/{a2 - Ci) and C[{x - 62)"^] = 
{Ni — A^o)/(^2 — Ci). If a2 7^ 62 then a partial fraction decomposition yields 

C[{x-a2)-\x-b2)-'] = {a2-b2r' C[{x-a2)-' -{x-b2)-'] = {No- Ni)/[{a2-ci){b2-ci)]. 

On the other hand if 02 = 62 then £[5*1 (x)] = yields C[{x — 02)""'^ + (^2 — ci){x — 02)^^] = 0, 
which implies C[{x — 02)"^] = (^2 ^ Ci)'^{No — Ni). We continue this argument by induction 
on n using C[x''Sn{x)] =0, < /c < n, n > 2. Thus for each new n,n > 2 we use 
^[^^(x)] = C[xSn{x)] = to evaluate £[njLi(2^ ~ (^j+i)'^ UkLii^ — &j+i)~^], N = n, M < n, 
N <n, M = n and N = M = n. 

The next result gives a recursive definition of £[a;" S'„(a;)]. 
Corollary 3.6 Set C[x^''Sn{x)] = Nn- Then the Nn 's satisfy the three term recurrence relation 

(3.4) Nn - Nn^i + A„ Nn-2 = 0, u > 2. 

Proof. Multiply the recursion (3.2) by x"~^, apply C and make use of the orthogonality 
relation £[x'^S'„(x)] = 0, < k < n. The result is (3.4). 

It is because of (3.4) that we need two initial normalizations A'^o and A^^i. In Example 3.2, 
which comes later in this section, we will see that Nn = /ti /€2 ■ ■ ■ i^n+i with /ti = 1/(1 — k,2)- 
Now Kj = Aj/(1 — Kj+i), J > 2 follows from (3.4). Since the continued fraction K^2 (^) of 
the example converges we find that (3.4) can be realized in practice by having 



Ki = (1 - K2) \ 



(3.5) /., = KZ, (y) , J > 1, 
and 

(3.6) Nn = K1K2 ■■■ Kn+l- 

We now come to a representation theorem for i?//-fractions. The polynomials of the second 
kind associated with the recursion (3.1) are generated by 

(3.7) Qn{x) - {x - On) Qn-i{x) + A„ (x - a„) {x - 6„) Qn-2{x) = 0, n > 1, 
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Qo{x):=0,Qi{x):=l. 

The ratio Qn{z)/Pn{z) is the nth convergent (approximant) of the continued fraction 
1 A2 (2; - 02) (z - 62) A3 {z - as) {z - 63) 



(3.8) Rn{z) 



Z — C\ — Z — C2 — Z — Cz 



which we assume becomes a finite fraction for z = a^ or z = 6^ for any k> \. This leads to 
the following definition. 

Definition. A continued fraction of the type (3.8) will be referred to as an i?//-fraction 
provided that it terminates in the cases z = a„ and z = bn, n > 1 . 

We will henceforth assume that Rn{z) converges to a function which vanishes at infinity 
and whose singularity structure is given by at most a finite number of branch cuts {Fj-jf^ and 
a denumerable number of poles {zj}'^ so that, in the generalized sense previously explained 
in Section 2, the following representation holds. 

da(t) 



(3.9) Rjj{z)=f^, zeC\V. 
Jr z — t 



Also, as before, we assume that V is minimal so that 



N 

(3.10) p = I U r, I U{^^}" 



and that the set of interpolation points {a„, &n}$^2; ^^^ [HI; i^ disjoint from the set of singular 
points V, that is 

(3.11) {a„,64~2nP = 0. 

We will also assume the large z asymptotics 

(3.12) Rji{z)^Ki/z = z-^ f da{t), 
with 

(3.13) '^i = Y_Y-T 

a convergent continued fraction. 
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Theorem 3.7 Consider the three term recurrence relation (3.1) and assume that the repre- 
sentation (3.9) and the conditions (3.10)-(3.13) hold. Then the linear functional C of Theorem 
3.1 has the integral representation 



(3.14) C[f] = J^f{t)da{t). 

Furthermore the normalization constants of Corollary 3.2 are realized by (3.5) and (3.6). 
Proof. We set Nq = ^i. it remains to show that Ni = ^i — 1 and that 

(3.15) f t'' Snit) da{t) = 0, 0<k<n. 

This will follow from the large z asymptotics in (3.16) of the Lemma below. 

Lemma 3.8 Under the assumptions of Theorem 3.3 we have the integral representation 

(3-16) . ^... .. " ^!l. ...^.„w . = / '-^da{t), 0<k<n,zeC\V. 



'■3 

where 



Ain-iM(^-«,)(^-&.)]^ir^^-^ ^r ^-^ 



(min) / \ / \-t^(mm)/ \ ^(min) , 



X,{z - ai){z - h)X^T \z) = {z- c^)Xr \z) - Xr \z) 

(see the remark before Theorem 2.3). In (3.16) X^^^\z) denotes the minimal solution to the 
recurrence relation (3.1). 



Proof. Here again we invoke Pincherle theorem and establish the representation 

A,(^-o,)(^-6i)AlT"'W *'--*' 
The asymptotic relationship (3.12) gives 

(3.18) ^\ X> ^n,z. 

The three term recurrence relation then implies 

X;™")(^)/Xi'"™)(^) ^ (1 - 1/K,)Z = K,Z, 
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say, with K2 := (1 — l//«i) and 

Xlr'-\z)/xiZ"\z) ^ z{l - A„/«:„) = ^n+iz n>2. 
This estabhshes the large z asymptotics on the left of (3.16) as 

k v-imin)( \ 



X,[Uni{z-a,){z-b,)]X}-r{z) 



"ni 



with Nn = Ki fi;2 ■ ■ ■ Kn+i and the k's are given by (3.1). To establish the equality (3.16) we 
follow the same procedure as in Lemma 2.4. This means that the singularities of (3.16) and 
(3.17) are the same but (3.16) has the additional weight factor t^ Sn{t), k < n. 

Note that the large z asymptotics of (3.16) yields 



/ t'' Sn{t) da{t) ^ z^-"" Nn, 0<k<n. 



The choice k = n = 1 yields A^i = ki k,2 = ki — 1 but the remaining choices give the 
orthogonality relations /p t'' Sn(t) da(t) = 0, < k < n as required by Theorem 3.3. 

Example 3.1: Chebyshev Polynomials of Type Ru. This case with constant coeffi- 
cients is given in |Q with some misprints which are corrected here. Consider the recurrence 
relation (3.1) of the form 

(3.20) Xn+i{z) -{z + V^)Xn{z) + -{z-a){z-b)Xn-i{z) = 0, a,b>0. 
This difference equation has solutions 

(3.21) X^iz)- 



±..A_ i{v^±v^){v^±Vb)Y 



For z ^ (—00, 0] the minimal solution is therefore 

(3.22) Xr"H^)=( ^^"^^,^^'^y 
with the square root branch chosen so that 

\z + Vab- {^/a + Vb)^/z\ < \z + Vab + {^/a + Vb)^/z\. 
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From Pincherle's theorem we then have 

^3 23) 1 {z-a){z-h)/A {z-a){z-h)/A 



z-\-\/ab— z + \fab — z + \fab — {^/z + ^/a){^/z + \/h) 



2 /"O {^ + Vh)^/^dx 



n J-oo {a — x){b — x){z — x) 

Corresponding to the more general formula (3.16) we also have, for < rn, < n, 

. . P^(^)X^-")(^)(v/i- v/^)(v/i- y/fe) ^ 1 ro i^+Vb)V^PUx)Pn{x)dx 
^ ■ ' (;z-a)"+i(z-6)"+i nJ-oo (a - x)"+i(& - a;)"+i(^ - a;) ' 

where 

.3 25^ p,^_ [(v^+v^)(v^+v^)r^-[(v^-v^)(v^-v^)r^ 

is the polynomial solution (no longer monic) with initial values P„i = 0,Pq = 1 

From the large z asymptotics of (3.24) and (3.25) we easily obtain the orthogonality 

.oofiX 2 fO (v/^+V&)P,„(x)P»(x)v/^C?X _^,2n+l^^ , ^^x 

^^•^^^ ^i-oo (a-x)"+i(fe-a:)"+i " I'^+iJVn- 

If we now introduce the rational functions 

/Q ov\ D / \ P2n{x) P2n+l[x) 

3.27 P2n a; := -— — , P2n+i x := -— , n > 

(a — x)"(o — x)*^ (a — x)"(o — x)'^+^ 

then (3.24) (with m = n and 2; = a) and (3.26) (with n ^ m) translate into the rational 
orthogonality 

2 fO {^+Vb)y^^RUx)Rnix)dx _ 2-^"a-V^ 
^^ ^ Tri-oo (a-a:)2(6-x) (^^ + ^6) '"^'"- 

Example 3.2: 2P1 Functions. A contiguous relation for a hypergeometric function of 
type 2P1 in [||, (45), p. 104] tells us that the recurrence relation 
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z{z-l){n-l){n + a-l-b) /^_n 

(2n + a - 1 - 6) (2n + a - 3 - 6) ^^-^^^> 

has solutions 

f'i'in^ vWf ^ f^V T{n+l)T{n+l + a-b) p/ . ^.^ x 

and 

f2W X /-z-lA" r(n + l)r(n+l + a-6) , , , x 

and a polynomial solution 

(3.32) Xi^\z) = Pn{z- a, h) := J % ^F.i-n, b - a - n;b - n;l - z). 

The asymptotic behavior of X^\z) and X^\z) as n ^ oo can be easily calculated and it 
yields 

(3.33) Xi^\z) ^ (2/2)"n-('^+''-^)/^ \z\ < 1, 

(3.34) X^^\z) ^ (^i^)\("+^-i)/2, \z-l\< 1. 
Thus the minimal solution of (3.29) is given by 

" ' ' \ .Y,f (j), Kj > 1/2, 
Exploiting Pincherle's theorem we obtain 

[ - ^^+"'\^^'^^"\ Fi(a,6;l + a;z), 7^z < 1/2 

(3.36) i?77(z) = I 

i -^^^z-''2Fi{l-a,l-b-2-b-l-z),nz> 1/2 

where 

(3.37) i?,,(^) = [^ - ci + K^Lsl^ (^ - 1) K/{z - c„)}]-\ 



(n + a — 1) (n — 1) (n + a — 1 — 6) 



(2n + a - 6 - 1) ' (2n + a - 1 - 6) (2n + a - 3 - 6) ' 
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We now have an explicit example of an Rn type continued fraction which converges except 
when TZz = 1/2. To calculate the absolutely continuous measure which is now supported on 
TZz = 1/2 we use [|, 2.9(33), p. 107], 

(3.38) ^'^ + ''-^\ l-zf~i^F^ia,b;l + a;z)- ^'^;l''~^K ~%F^il-a,l-b;2-b;l-z) 
a (o — 1) 



z-"" il-z 



,^{l + a-h)T{a)T{l-h) 



V{l + a-h) 

Thus we have established 

il + a-h)Via)Vil-h) fi+^'^t-'^il-tf-^ , , , 

3.39 Rn{z)= ^ ^ .'^^^ \. ^ I ^ ^rft,7^^^1/2, 

27rzr(l + a — o) Ji-joo z — t 

with the correct 0(1/2;) asymptotics of the present section when TZ{a — b) > which is 
2.1.4(17), p. 63] 

(3.40) Rii(z) ^ r~^tM z-\ n(a - 6) > 0. 

Y a — b J 

This justifies the assumption (3.12) since we may rewrite the A„'s in the form 

, , 1 l-(a-6)2 

(3.41) A„ = -+ ^ ^ 



4 4(n + (a-6-3)/2))(n+(a-6-l)/2) 
and a result in [0 yields 

1 Ao Ai CL — b + 1 „ , 

3.42 K^ = - ^ ^ = — , 7^ a - 6 > 0. 

1— 1 — 1— ■■■ a — b 

We now derive a biorthogonality relation for certain rational functions using the orthogonality 
relation (3.15) with the special case fc < n, 2; = of the integral formula (3.16). We set 

PJx; a,b) 

(3.43) Unix; a, b) := -^ — 4^;-, Vn{x] a.b) := f/„(l - x; -b, -a), 

[x — 1)" 

(3.44) g{x; a, b) := x'^-^ (1 - xf-^ /2m. 

The orthogonality relation (3.15) now takes the form 

(3.45) / t-"'Unit;a,b)g{t;a,b)dt = 0, < m < n. 

21 



Taking the complex conjugate of (3.45), replacing a by —h and h by —a, taking into account 
that t = 1 —t and using (3.16) (with z^ replaced by Pk{^ — z\ —b, —a), k = n and 2; = ) we 
see that 

provided that a, 6 7^ 0, a 7^ —1, 6 7^ 1, 7^(a — &) > 0. 

The biorthogonal rational functions in this example differ from the ones considered by 
Askey in . Both systems of rational functions are biorthogonal with respect to a Cauchy 
beta integral weight function. 

4. The 403 Functions. In this section we provide two examples of biorthogonal rational 
functions which have explicit representation as 403's. The 

three term recurrence relations and the continued fractions associated with these systems 
come from contiguous relations for 4^3 and §07 function. The 4^3 contiguous relation is the 
same one used by Askey and Wilson [^] and that led to the Askey- Wilson polynomials. The 
difference between our examples and the Askey- Wilson polynomials is a reparametrization 
which converts the three term recurrence to one of type i?//. 

The monic Askey- Wilson polynomials are defined by 



4.1 Pn{x;a,/3,j,6\q):= ^ 403 . . 

{2a)"- {al3-fd/q;q)2n \ aP,a-f,a6 

where 

(4.2) x:= -{u + l/u). 

The PnS satisfy the three term recurrence relation 

(4.3) Pn{x) -{X- an-l) Pn-l{x) + A„_2 B^-l Pn-2{x) = 0, 



an = -{a + 1/a) -An -Br, 



_ (1 - a/37(5g"-i)(l - a/5g")(l - 079") (1 - aSq") 

An '■ — 



Q^Q 



'-n 



2a (1 - a/575g2n-i) (^ _ af3'^6q^") 
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_ a(l - g")(l - /37g"-^)(l - /5(5g"-^)(l - 7(5g"-^) 
"~ 2 (1 - a/37(5g2„-2) (i _ a/5^5g2„-i) ■ 

In [jTB[ the minimal solution to the recurrence relation (4.3) for n > has been shown to 
be 

(4.4) Fi--)(x) = | ^"^''^' ''''>^' 
\ i^n(l/«), |«| < 1, 

and the F„'s are the functions 

(2m)-" (a/37(5g2"-i, ag"+7u, /3g"+VM, 7g"+VM, 5g"+Vu; g)oo 



(4.5) FM ■■- 



(g"+^, g"+2/M2, a/3g", a7g", a(5g", /37g"', /3(5g", 7(5g"; q)c 



\ aw pu 7M OM / 

From (4.4) one can calculate, using Pincherle's theorem, the Askey- Wilson weight function 

(A 6) w(x) ■= "*" (m^ 1/m^ a/3, ^7, a6, p'y, (36, -f6, g; g)oo 

27rVl - a;^ (aw, a/u, Pu, (3/u, ju, j/u, 6u, 6/u, aP'j6; g)oo ' 

normalized by 

1 

w{x) dx = 1, |a| < 1, \j3\ < 1, I7I < 1, \6\ < 1. 

The system of biorthogonal rational functions recently given by Al-Salam and Ismail []1| can 
be related to the above Askey- Wilson case through a change in parameterization. This will 
give us our first explicit 4^3 example of an i?//-fraction. 

Example 4.1: Biorthogonality on the Unit Circle. Let us use the new parameters 

a, b, ti, ^2 and the variable z given by the replacements 

(4.7) a = q-^/%^, (3 = i'^a^z, 7 = hq^'^j^z, 5 = q-^1%/^, u = q'^'^^z. 
Then after renormalization (4.3) becomes a recurrence relation of type Ru, namely 

(4.8) Pn{z) -{Z- C„) Pn-l{z) + \n{z - a,) {z - 6„) Pn^2{z) = 
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with 

(4.9) A„,+i 



at2g"(g" - 1)(1 - a6tit2g"~^)(l - 6^2?""^) (1 - tit2g"-2)(l - abq"){l - atiq""-^] 
4 ^(1 - a6tit2g2«-3) (1 - a6tit2g2"-2)2 (1 - ahtit^q^'^-^) u^ u„+i 



_i = 6tig"-3/2, 6„^, = g-"+V2/(at2), c„ = -t;„K, 



-1/4 



Mn+1 



2(1 - abtit2q'^'''^)il - abtitiq 



2n) 



X {1 - g""i[(l + ^^^q'-)(qt, + ag^ + at.h + abhq) ""''''' 



q 



abtit2 „,^ ,, q g^,-,, 

^i^g"(l + g)(t2 + ga + H f)]} 



,1/4 



Wn+1 



2(1 - a6tit2g2""^)(l - a6tit2g 



2n\ 



X {1 - 5»-"|(l + ^^^,2")(6, + ii + ait, + 



q q q t2 a 



Note that, in order to calculate the middle coefficient in (4.8) and arrive at values for c„, Un 
and Vni it is very useful to use the form of the Askey- Wilson three term recurrence relation 
given in [^, (1.24), (1.27)]. This is true also for the next example. 

The polynomials {P„(z)} are now given by 



(4.10) Pn{z) 



g"/^(at2zg^/^ 6^2, tit2/g, abtit2/q] q)n 
{2t2)''{abtit2/q; q)2n Uk=iUk 



X 403 



g-", abtit2q''-\ t2zq-^/^, ta 
at2zq^/'^, bt2, tit2/q 



q,q 



From (4.4) with a renormalization factor z"'^'^/ YVk=_i Uk we obtain the minimal solution to 

(4.8) 



(4.11) X^™")(2) 



XW^z), |.|>|g|V2 
X(2)(^), |^|<|g|V2, 
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with 






/ n+3/2 3/2 1/2 2 \ " 

\ ^ t2Z UZ t\ 



k=-l 



and 



(4 13) X^'Hz) - {q-'"z/2r{ahUhq^^-\ hzq^^^'^ azq^+^l^ bq-+\ tig"; g)oo 

" (g"+^zg"+3/2^Q^2zg'^+V2,6t2g",tlt2g"-^a6g"+^atlg",6tlg"-V2/;^;g)<. 



V t2 a b t, J ^i\ 

With the above minimal solution, equations (3.9) and (3.17) imply 

(,.,. p ( . 2mq-'/\l-ahht,/q){l-qy^z) 

^' ' ''^' {l-b){l-h/q){l-t2zq-y^){l-azqy^) 



Cn 1 a '7 fi ' 7 \ 

^g'/'; g, f, -, ^, ^; a6tit2/g , /or |z| < \q\^i'' 
t2 a b ti I 



(4.15) Rn(.) 2u,,"Hl-abt,tM(l-,^'y. 



z{l - t2)(l - aq){l - 6gVV^)(l " hq-^/yz) 

/„3/2 3/2 1/2 2 \ 

X sW^r — ;g,^,— ,|,f;a6tit2/g , /or |^| > Igl^/^. 

\ Z T2Z OjZ t\ I 

An application of Cauchy's theorem gives us the integral representation 

a'{t) 



(4.16) Rn{z) = f 

J\t\ = \a\i 



dt, 

|tl = |g|l/2 Z-t 



1/2 



with a'(t) as ^ times the difference of the boundary values of i?//(z) in (4.11) as \z\ -^ \q\ 
The identity [^, (III. 37), p. 246] then gives a'{t) exphcitly. The result is that a'{t) is ui/t^/^ 
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times the Askey-Wilson weight (4.6) but with the new parameterization given in (4.2) and 
(4.7), with z replaced by t. In order to isolate the sjTiimetric terms in a' we chose to write it 
in the form 

4.17 a{t) = — ^- rr/t 

where /(t) is 

{q^^H, q^''^/t, at2q^'H, btiq^^'^/t, bt2, ati, abq, tit2/q, q; q)oo 



(4.18) fit) 



{aq^/H,bq^/^/t,aq,bq,ti,t2,q ^/%t,q ^/Hi/t,abtit2\q)oo' 
Note that /(t) is symmetric under the transformation 

(4.19) {a,bMM,t) -^ {b,a,t2,h,l/t). 

We now derive the biorthogonality relation of Al-Salam and Ismail [0]. In the case under 
consideration our orthogonality relation (3.15) becomes 

(420) ^^/ t^p^{t){i-bt,q-yyt)mdt ^oo<^<^ 

^- ' nqy'J\t\=My-t-ibhq-yyt;q)^UU[t-q-'^'^yiat2)] ' " 

We follow the notation in |l[] and define rational functions r„ and s„ by 

/ , X / q-'',abtit2q^-^,q-^'H2Z,t2 

(4.21) r„(2;) =r„(z;a,6,ti,t2) :=403 . 1/2 r,. . . / ^' ^ 

y at2q ' z,bt2,tit2/q 

and 

(4.22) Sn{z) = Sniz; a, b, ^1,^2) := r^iz; b, a, ^2, h). 

There is a slight difference between the above definition of s„ and the definition of s„ in |jl[] 
due to complex conjugation. Now the relationship between P„ and r^ is 

(4.23) r^{z;aAti,t2) ^"^"^^^ 



{at2q^/'^z;q)n' 
where the normalization constant Cn is given by 

{2t2q-'/T {abt,t2/q;q)2nUl=iUk 



^n 



{bt2,tit2/q,abtit2/q]q)n 
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We may now rewrite (4.20) as 

(4.24) / ^ r„,(t) Sk{l/t) fit) - = 0, k<n. 

If we deform the contour in (4.23) to the unit circle |t| = 1 and combine this new integral 
with its complex conjugate (for which t = 1/t) then perform the replacements (a, b, ti,t2,q) -^ 
(b, a, ^2, ti,q) while taking into account the symmetry of / under the transformation (4.19) we 
obtain the biorthogonality relation 

f dt 

(4.25) / rfc(t;a, 6,ti,t2)sn(lA;a,&,ii,i2)— = 0, k ^ n, n,k>0. 
■J\t\=i t 

The value of this last integral when n = fc is deduced from a modification of (3.16) with k = n 
and z = fetig"^-*^/^ < 1. Namely, 



(4.26) 



P„(l/;2)X(™")(z) f rP„(l/t)P„(t) a'{t) dt 



Ai[n"i'(^ - «,)(^ - &,)]^ii (^) ^1*1=1^1^^^ i^l=2it - «.)(i - b,)] (z - t) 



where P„(-2) = P„(x; 6, a, ^2,^1)- We shift the above contour to |t| = 1 and put z = btiq" ^/^. 
At this value of z the sWj in XJ^'^"''{z) is summable 

using ip, III. 23, p. 243]. After a lengthy calculation we arrive at 
(4.27) ^f r^{t)sMt)m'^'- it^t^/^ri<l,abq,abt,t,q-~^;q)^ 



2tt J\t\=i t {tit2/q;q)n{abtit2;q)2n 

in agreement with |^. 

The condition for the absence of the discrete spectrum and the validity of the integrals on 
the contour |t| = |g|^/^ is |a|, |6|, |ti|, |t2| < 1- The condition for shifting the contour to |t| = 1 
is to have \a\, \b\ < \q\^^^'^ and |ti|, |t2| < \q\^^'^- Thus our overall condition for deriving the 
biorthogonality is \a\, \b\ < 1 and |ti|, |t2| < Iq]^^"^- 

Example 4.2: Biorthogonality on the Line. We next use a different set of parameters 
^1,^2,^35^4 and variable z. We set 



(4.28) a = VhUh/q, (3 = -qe^/Vtsk, 7 = qe'^/Vtsh, S = t2g"Vv^ 



4, 



u = —yt^/t^, z = sinh^. 
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The explicit representation (4.1), after renormalization by 11^=1 ""A; is now 

(-tie^, tie"^, tit2htiq-^, -titsg"^ q)n 



(4.29) Pn{z;tiMMM\q) 



(2tiv/t^/g)" (-tit2/g2;g)2„nLiWfe 



X403 






9,g , 



where «„ is given below in (4.33). 

The three term recurrence relation (4.3), after renormalization is now of type Rn and is 

(4.30) Pn{z) -{Z- Cn) Pn-l{z) + Xn {z - a„) {z - 6„) P„_2(^) = 



with 

(4.31 



) A, 



n+l 



tit2g'"-'(l - g")(l + tit2g"-3)(l + g"+V^3t4)(l - hhhUq'^-*) 



;i + tit2g2"-2) (1 + tihq^^-^y (1 + tit2g2"-4) 



UnUn+l 



(4.32) a„H.i = i(t2g"-'-g'-"A2), 6„+i = i(tig"-i-g-"-Vti), 



Cn+l 



-^^n+1 + \{\Jh/U + yW^) 



Mn+1 



where «„ and f „ are given 



(4.33) M„+i = [(i-tit2g'""')(g' + tit2t3t4) + g"(i + g)tit2(g-t3t4) 



;i + M2g2«-3)(i + tit2g2"-i; 



(4.34) w 



n+l 



-[(1 - tit2g'"-^)(t3t4 - q') + g"-'(l + q)ihhhU + g') 



2n-l^ 



(4.35) 



2v/titl(l + tit2g2"-3(i + tit2g' 

In the case under consideration the continued fraction is the Ru fraction 
1 A2(-2- a2)('2- ^2) 



Z — Ci 



Z - C2 
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When the continued fraction converges it will converge to F{z), where 

2mi il-qU/t3)il + ht2/q^) 



(4.36) F{z) 



^ts/U (9^4/^3, -titi/q, qe^ts, -qe V^3, -^2^4/^^)00 
X W{qU/ts; q, -q^/ht^, Ue'^ , -Ue^ , -q^ /t2h; tits/g^) 



with 

crq a q 

W{a- b, c, d, e, f; t—t^) ■= (&, c, d, e, f; g)oo8W^7(a; b, c, d, e, /; -r-r^) 
bcdej ocdej 

The singularities of (4.36) are at z = z„, 

Zn = l{t3q-''-'-q''+'/h), n = 0,l,---. 

given by the zeros of {qe^/t^, —qe^^/ts; g)oo- The residues at these points may be calculated 
explicitly because the contribution from W is in terms of a very well-poised 605, which is 
summable. After a straightforward but rather lengthy calculation, we arrive at the Mittag- 
Leffler expansion 

00 
(4.37) F{z) = J:-^ 
fc=0 ^ ^^ 

with 

ui^/t^A {tit2hUq~^, qh/ts, ^2/^3, 9^4/^3; g)oo 



(4.38) LOk 



qik+i (-t^t2/q, -hU/q, -t2U/q\ -q^/tl q)oo 
,, (-g /tih,-q /ht3,-q /ht4,,-q /%;q)k i^ . ■2k+2 i+2\(. + + + \k 

r~n7 — 4~77 — 4~Ta — r^; u + g /t3J(^tit2t3t4j • 



iqti/t3,t2/h,qt4/ts,q;q)k 
The resulting orthogonality (3.15) is now given by 

A;=o lli=il^fc - aj+ij(^2;fc - Oj,+ij 
We now derive a rational biorthogonality relation. Let us include the factor 

Zk-a2 = ^jq-'-'il - t2q'/t3){l + q'^Vt2ts) 
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in the weight to obtain the new weight 



Tk 



UJk 



Zk — «2 

2t2Uiy/hU jtititshq'^ , gti/tg, gta/tg, gt^/t^j-, q)^ 
q^^iq^ + t2h) (-tit2/g, -hU/q\ -t2U/q\ -qV^l Q) 
i~q^/tit3,-q'^/t2t3,-q^/ht4,-q'^/tlq)k , 2fc+2 /^2 



X 



{qh/ts, qt2lh, qti/h, q; q)k 



l + q''^'lti){ht2ht,f. 



Note that this has k dependent terms which are now symmetric in ti and t2- This symmetry 
together with (4.39) yields the biorthogonahty relation 

oo 

(4.40) Y.^^^^k)Rn{zk)rk = ^, m^n, 

k=0 

where 



(4.41) Rn{z) 



Rn\Z] ti,t2,t3,t4j 

-he^, he'^, tit2hUq~'^ 



4(A'3 



q^ 



Rn\z) — Rn{z', t2, tl, ^3, ^4 



and from (4.29) we have 



(4.42) — 



Pn{z) 



-irq 



n n(n+l)/2 



(t3t4 



-n/2[ 



-tih/q^.tihhtiq ^;q)n 



U]=i{z~b,+^) ' '^ '■"-' {-ht2/q';q)2nm=iUk 

For the case m = n we use (3.16) with k = n and z^ replaced by 
Pn{z;t2,ti,t3,t4^) and z = an+2- This requires the evaluations 

q-<^-^)l\-t2h/q,-t2U/q-q)^ 



RJz) 



(4.43) Pnian+2;t2,ti,ts,t4 

and from |, (111.24)] 



-2t2Vhk)^ 111=1 uk 



n+l 2 



(AAA) 8W^7(g"+'t4A3;g"+\-g7M3,t4e~S-t4eS-g7t2t3;tit2/g)|e«=t2g-i 



? _«3 



(g"+'t4/t3, -kU/q, g"t2/t3, -g'"-^tit2; g) 
(g2"+it2/t3, -g"tit4, gt4/t3, -tit2g"-2; g)c 
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The final result after another lengthy calculation, is 

(AA^\ V R r. IR (. ^ (-gV^i^3, -gVt2t3, -qVhU, -q^/tl q), 1 + g^^+Vt^ , 

(4.45) Ei?.(z.)i?„(-.) (gt,/t3,gt,/t3,gt,/t3,g;g), TT^AT^'^'^'^'^/' ^ 



(W3Wg^^"^ + '^'^/^' (V/t3t4,g;g), 



1 + titag'"-' {-tit2lq\ ht2t3U/q^; q)n 



(-tita/g, -tiU/q, -t2U/q, -q^/tl; g)oo 

lOD 



{qh/h, qh/h, qti/ts, tit2hti/q^] q). 



with 



Zk = 2 (^sg"''"' - g'^+'As), A; = 0, 1, ■ ■ • , \tit2ht4\ < g'. 

Another biorthogonality relation can be obtained by interchanging ^3 and ^4, which corre- 



sponds to taking u = —Jt^/ts. Both of these are special cases of the biorthogonality relations 
derived in [|l^, see (1.21), (3.14) and (3.15) in [l^, with a = q/t-^ or a = q/t^^. However, if 



in (4.28) we had used (3 = qe^ / y/t^t^ and z = cosh^, then we would have arrived at a similar 
biorthogonality with mass points at 

^fc = ^(^3g-'-' + g'+'A3),fc = o,i,---. 

We do not know where this fits into the scheme of things other than to suggest that, with 
discrete orthogonalities at this level, there are always two families, one associated with z = 
sinh^, and another with z = cosh^. 

5. Biorthogonality on [—1, 1]. We present two final examples by yet another modifica- 
tion of the Askey- Wilson recurrence relation. Here an i?//-type continued fraction is used to 
obtain two g-beta integrals, one which is new and one which is in 



2% . However the biorthogonalities obtained are not given by the continued fraction 
methods of the previous examples. Instead we use the "attachment" method explained in [^ 
which goes back to Andrews and Askey and was used in |]5|, [|I|]. 

In (4.3) we apply the parameter replacements 

(5.1) f3^f3/u, i^(3u, 

31 



with z = {u + l/u)/2. The recurrence relation again becomes an i?//-type recurrence relation 
(4.30) but now with 



a/3=^(5g^"-^(l - g")(l - a/j^^g"-^)(l - a^g"-^)(l - /j^g"'^) 



(5.3) a„+i = ^(/3(5g"-i + gi-7/3<5), &n+i = ^iaPq''-' + g'^V^/^), 

(5.4) Cn+l = Vn+l/Un+1, 



where «„ and Vn are given by 

„_i (1 + aP^6q^''-'){q + a6) - g"-i(l + g)(g + /32)a(5 



(5.5) Un+l = 1 - 



(5.6) t;„+i = (a + (5)g' 



(l-a/325g2n-2)(l_^^25^2„) 
,_1 (1 + a/52(5g2-l)(g + /52) - g"-l(l + g)/32(g + a^) 



2(1 - a/325g2«-2)i _ a/52,5g2") 
The new minimal solution of the recurrence relation is given by 

(5.7) xi-\z) = { f-^f' i;i>\ 

Gn{l/u), \u\ < 1, 



n 

x,W^{q-'-^/u^-q-+\q/au,q/(3,qu-^/(3,q/5u-a[3Hq'^-^)/ \{ u^. 

k=-l 

An application of Pincherle's theorem now yields the Stieltjes transform 
^ f{x)dx 2 {l-q/u^){l-af3^6/q) 



(5.9) 



1 z-x u{l-a/u){l-p/u'^){l-p){l-5/u) 

X sW7{q/u'^; q, q/au, q/(3, q/(3u'^, q/6u; ap'^6/q), 



with z = |(-u + I/m), !«! > 1, and 

1 (e2^^e-2^'^,a/3e^^a/5e-*^/3(5e*^/5(5e-^^a(5,/3^g;g)oo 1 



(5.10) f{cose):-- 



27r (ae^^ae-^^/5e2^^/5e-2^e,(5e*^(5e-^^/3,/3, a/525; g)oo sin^' 
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The integral (5.9) seems to be new. 

The large z asymptotics of the transform formula (5.9) produces the following result 

(5.11) ^£/(cose) sin^ci^ = ^^^20i(g,g//3;g/3;g,a/5'Vg), 

valid for max{|a|, |/?|, |(5|, |Q;/?^(5/g|} < 1. 

We next explain where (5.11) comes from in the theory of orthogonal polynomials. Recall 
that the continuous g-ultraspherical polynomials {C„(a;; /3|g)} have the generating function 

(5.12) t Cn(co. 0; P\,)t" = "|g!;fl,2)^ - 

and satisfy the orthogonality relation 

1 PTT (p2ie -2ie. \ 

(5.13) — / CUcos9;P\q)Cn{cos9;P\q) )^ ^^; J^°^ d9 

{(3,q(3;qU{(3';qUl-(3) 

(g,/5^;g)oo(g;g)n(i-/3g") '"'"' 

0, 0. It is easy to see that (5.12) and (5.13) show that the left-hand side of (5.11) is 

= E ^/ cUcose-,mCn{cose-,f3\q) .X^,,;'^J, dO 

^7lo 2vr Jo (/3e^*^/3e ^*^;g)oo 

f^. ,.n (/3,g/3;g)oo(/3^;g)n(l-/3) 

{(3,q(3;q)^ ^ 2 



q,P^]q)c 



20i(/? ,(3;q(3;q,a6). 



Now the Heine transformation p, (III. 3)] reduces the extreme right-hand side above to the 
right-hand side of (5.11). 

Example 5.1 Chebyshev Rational Functions. The relationship (5.11) gives an integral 
representation for a basic hypergeometric function of the type 20i- In this generality the 
integral representation (5.11) does not seem to lead to orthogonal or biorthogonal functions. 
The special case f3 = q degenerates into the elementary integral result (1.3). To find rational 
functions biorthogonal with respect to the integrand in (1.3) we now use the attachment 
method. 
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Consider functions of the type 
(5.14) 9.(oosO;a.^) := g ^__j-_-5-^a,,,. 
where an,k are to be determined. Set 

^ ^ "'^■~ Trio ^"^""^ '"' \g<5e^^g5e-^;g),■(l-2acos^ + a2)(l_25cos^ + 52) 



sin^ 9d9 



Thus 



"■n,] 






Vl — x"^ dx 



fc=o (^' ^)fc ' ^ ^^ (^ ~ 2aq'^x + Q;^g^'^)(l — 2(5g-5'x + 5'^q^^) 



1 y^ i^g ,aoq-';q)k ^ 
1 - a(5g^' ^ (g, a6q^+^; q)k '^'^' 



Recall that the g-analogue of the Pfaff-Saalschiitz theorem p, 11.12)] is 



(5.16) 302 , 



q "',a,b 
c, abq^~"'/c 



Q,Q 



{c/a, c/b; q)n 
(c, c/ab; q)n 



We would like to choose an,k to make /„j vanish for < j < n. In view of (5.16) we let 
On.fc = q^i, aSq"", q)k/{(yS; q)k and we obtain 

iq~'',q~^;q)n 



"'^ {l-a6q^){a6,q-''-^/a6)n' 
It is now clear that In,j = ii < j < n. After some simplification we find 



'"J 



[(_a5)"g"("-i)/2(g; q)l/(a6- qhn+i]hn- 



Thus the gnS of (5.14) are given by 

q~^ (ySq^ cxc (y.c~ 
(5.17) ^n(cos^;a,5) =403 I ' '. ' _,. 

ad, gae , gae 



g,g 



and satisfy the orthogonality relation 

('•'') I L ^-^"' ^' "^^"^"' "' ^^ (1 - 2ax /j) (1 -^2^ + ^^) 
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^ (a(^)"(g,g;g)n ^ 

We find the biorthogonality relation (5.18) very surprising since the 4^3 function in (5.17) 
is not balanced. Recall that a basic hypergeometric function (1.9) is balanced if r = s + 1 and 
qaia2 ■ ■ ■ as+i = 6162 ■ ■ -bg- Only balanced 403's with argument q satisfy three term recurrence 
relations. This concludes this example. 

In the previous examples of pairs of biorthogonal functions the second family was obtained 
from the first by symmetry and permutation of parameters. In the next example the members 
of the pair of biorthogonal rational functions are not related in this manner. 

Example 5.2 A Nonsymmetric Case. Following previous examples we now consider 
the g-beta integral obtained by absorbing the factor x — 02 into the weight function. This 
means evaluating (5.9) at z = 02, that is u = l/ajS, and this is exactly the case that makes 
the right-hand side of (5.9) summable. This gives the result 

J_ r (e^'^ e-^-^ a(3qe\ a(3qe-^\ (35^', (35e-^\ aS, (3\ q; qU ,, 
^ ■ ' 271 Jo (ae^^ae-^^/?e2*^/?e-2^^5e^^(5e-*^/?,g/?,a/?2^;g)oo 

max{|a|, |/3|, \6\} < 1. 



l-a^(3 



After we showed an earlier version of this paper to Mizan Rahman he pointed out that 
this g-beta integral is not new. It was obtained previously by him in his interesting work |^ 
using a different method. It may be of interest to mention here that (5.19) is also the special 
case 7 = g/3 of 

J_ r je'"', e-^-^ a^e^', a^e-^', /3fc-^ /3fc--^ a6, /j^ g; q)^ ,. 
^ ■ ^ 271 Jo {ae(>,ae-^o,pe^i^,(3e-^i^,6e^,5e-i^,p,qP,ap^;q)^ 



{(3,j,a'^'j,(3^a6;q)oo , / a^P,a6,qp/j 



3(P2 



{q,l3'^,a'^f3,a6;q)oc \ a^'y,P^a5 



Q,l 



for max{|a|, |/3|, I7I, \S\} < 1, which will appear elsewhere. In view of (5.12) and (5.13) for- 
mula (5.20) is equivalent to the evaluation of the connection coefficients for the continuous 
g-ultraspherical polynomials. L.J. Rogers solved this connection coefficient problem in 1893, 
seei. 

One can use the attachment method of the previous example to derive a system biorthog- 



onal with respect to the integrand in (5.19). Rahman also has done this in |23]. For com- 
pleteness we repeat the calculation here. 
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The g-beta integral we are working with is 
(5.21) 
where 



w{x; a, (3, 6)dx = — , 

-1 {l-a^(3){a6,(3\q;q)oo 



(5.22) w{cose;a,p,6):-- 



Choose 



{e^'\ e-^'\ qa(3e\ qa(3e-'\ (36e'\ (36e-'^; g)oo 
2tt sin e {ae\ ae-^\ ^e^^ , (3e-^^<^, 6e^^, 6e-^^; q). 



(5.23) ^n(cosg;a,A^) = E ^^^^^^,^^^^_,.^^^ a„,.. 

Using the g-analogue of the Pfaff-Saalschiitz theorem, (5.16) we put 



(ln,k — q 



{a5]q)k 



and find 

(5.24) 



(cxc cxc q) ■ 
^(-os 9; a, P, 5)^„(cos 9; a, P, 5) ^^^^^^^^ .^^^e"^^; g), 



sin 9d9. 



iP,qf3,aP'q-6;q), 



(/3^g;g)r 



;i - a2/3g2n) (^^^n^ ^2^ ^. ^)^ (a/^s^g^, gi-«/a5; ql 



-Sj,n, 3 < n. 



Let 



(5.25) V^n(cos 9; a, [3, 6) := 4^3 



and 



g-",a/32(5g"-^(5e^^5e-*'' 



a5, (35e'\ p5e 



-id 



Q,q 



n ( „—n 



ae*", ae *"; q) 



(5.26) 0„(cos ^; a, P,5):= g (^^^^^.,^ ^^^^_.,. ^)^. 



^Q^h 



q'On,j. 



Thus we have proved that the rational functions {0„} and {ipn} satisfy the biorthogonality 
relation 



(5.27) / (f)m{x;a, P,S)^n{x;oi, P,S)w{x;a, P,5)dx 



{(3,q(3,a(3^5q'^;q)^ {(3^,q;q)n 
{a5,l3\q-q)^ (1 - a^/^g") 



{-a5Y bn,n Sm,n, fU < fl. 
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We now select bn,j in order to have full biorthogonality, that is to make (5.27) hold for all m 
and n. In order to do so we need to compute the integrals Jn,k, 

Jn,k-= / (pnicosO; a, (3,6) ' — '--^-w{cose;a, f3,S)smede, 

for k < n. Substituting for 0^ from (5.26) and using the g-beta integral (5.21) we get 
^ {(3, g/3, af3^6q''; q)^^ " (g"", a^g^ q)j q^ b^j 

In order to apply (5.16) we must choose bnj as 

{l-a^(3q^^){a(3'^6q''-^;q)j 



(5.28) bn,j 
hence 

'Jn,k 



[1 - a^p) {a6; q)^ 



{P,qP,aP'5q'^;qU (g-^ g^'V^ ; ?) 



(/52, aSq''; q)^ {a6, q^-^-'' /a[3'^5] q)oc ' 
Thus (5.26) becomes 

('5 29) d) (cos^ a B S)= 6 ( ^~"' «/^^'^^"~^ "^'^ «e~'^ ^«v^' -^«v^ 

I a6, qal3e^^ , qa(3e~^^ , a\f]3, —ay/P 

This establishes the biorthogonality relation 



(5.30) / (prni^ia, P,5)ijjn{x;a, P,5)w{x;a, P,5)dx 



q,Q 



{13, qP, aP'Sq^; q)^ {(3^, q, ; g)„ (1 - a^/J^g^-i) (a6) 



■'m,ni 



(1 - a^f3){a6, (3\ q; q)^ (1 - a^p6q^-^){a6; g)„ 

valid for |a|, |/3|, \5\ < 1. 

Special cases of the -j/^n's have appeared in the literature as special cases of the Askey- 
Wilson polynomials. The case /3 = gives the Al-Salam-Chihara polynomials, [Q, [^, as can 
be seen by comparing (5.26) and (5.29) with (4.1). The case a = 5 = of (5.21), (5.22), (5.26) 
and (5.29) gives the g-ultraspherical polynomials of L. J. Rogers 0, but the justification 
is not as straightforward. It can be justified however by using the generating function H^ 



Y {q,cd;q)n \ ab, ac, ad 



Q,Q 
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49 



\he 



ah 



1,(1 I 201 



ce , de~ 



cd 



Q^q 



In a private communication R. Askey mentioned that he can let a ^ in (4.1) and show 
directly that \ima^oa~'^Pn(x;a, P,'y,6\q) exists and find its value. We do not know how 
to compute the limit \im.a->Qa~"'(f)n{x;a,P,6\q) directly but we suspect that Askey's direct 



argument may work here. In [^ M. Rahman made the transition to the a = 6 = case by 
using a 302 representation for the g-ultraspherical polynomials, writing 0„ as a combination of 
403's, applying Sear's transformation for 403's to the ipn and the 403's in 0„, and then taking 
the limit a,6 ^ 0. 
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